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result from the long-range order in the various phases are determined, and their consequences for 
electronic properties, in particular the specific heat, the single-particle relaxation time, and the 
electrical and thermal conductivities, are derived. Various aspects of the ferromagnetic limit, and 
qualitative differences between the transport properties of helimagnets and ferromagnets, are also 
discussed. 
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A. Dzyaloshinsky-Moriya magnets 

Helical magnets are systems in which long-range mag- 
netic order takes the form of a helix or spiral, such that 
in any given plane perpendicular to a preferred direc- 
tion there is ferromagnetic order, but the direction of 
the magnetization rotates as one goes along the pre- 
ferred axis. The pitch vector q of the helix points in 
the preferred direction, and its modulus q = \q\ is the 
pitch wave number, with 2jr/q equal to the helical wave- 
length. One mechanism for stabilizing this type of order 
over the homogeneous ferromagnet was pointed out by 
Dzyaloshinksky 1 and Moriya^ It relies on the spin-orbit 
interaction, which can lead, for certain lattice structures, 
to a term of the form M ■ (V x M) in the Hamiltonian or 
action, with M the magnetic order parameter. The pres- 
ence of such a chiral term implies that a homogeneously 
magnetized state can always gain energy by means of a 
nonzero curl of the magnetization, which leads to a heli- 
cal ground state. In a rotationally invariant system the 
direction of the pitch vector would be arbitrary (analo- 
gous to the arbitrary direction of the magnetization in 
an isotropic ferromagnet). However, the underlying lat- 
tice leads, via crystal-field effects, to a pinning of the 
helix in certain directions determined by the lattice. The 
crystal-field effects are due to the spin-orbit interaction, 
as is the Dzyaloshinsky-Moriya (DM) mechanism for he- 
lical order itself. Since the spin-orbit interaction is weak 
on an atomic or microscopic scale, with a dimensionless 
coupling constant <? so <C 1 , this leads to a hierarchy of en- 
ergy or length scales that can be classified according to 
their dependence on powers of <? so . Finally, an external 
magnetic field provides an additional energy scale, and 
couples to the helix via an induced homogeneous com- 
ponent of the magnetization. This leads to a rich phase 
diagram that is the topic of the present paper. 

1. The phase diagram of MnSi 

A very well-studied helimagnet of DM-type is the 
metallic system MnSi, which we will concentrate on; 
a very similar, but less extensively studied material is 
FeGe. 3 In zero magnetic field, and at ambient pressure. 
MnSi displays long-ranged helical order with 2ir/q m 
180 A below a temperature T c 30 K. 4 MnSi crystallizes 
in the cubic B20 structure, and the helix is observed to 
be pinned in the (1, 1, Indirections.^ The corresponding 
space group is P2i3 X In a magnetic field H, the helix is 
superimposed by a homogeneous component of the mag- 
netization, which leads to the so-called conical phase.— A 
magnetic field in the [0, 0, Indirection tilts the helix away 
from (1, 1, 1) until the q-vector aligns with the magnetic 
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FIG. 1: Schematic phase diagram of MnSi in the H-T plane 
showing the helical, conical, and A-phases, as well as the field- 
polarized ferromagnetic (FM) and paramagnetic (PM) states. 
See the text for a description of these phases. 



field at a critical field strength H c \. With increasing H 
the amplitude of the helix decreases, and finally vanishes 
at a critical field H C 2, where the system enters a field- 
polarized ferromagnetic phase. Inside the conical phase 
at intermediate fields near T c there is a region known 
as the A-phase4 This was thought to represent a he- 
lix whose pitch vector is perpendicular to the magnetic 
field^ii^ but recently has been interpreted as a topolog- 
ical phase where three helices with co-planar q-vectors 
form a skyrmion-like structure.— The schematic phase 
diagram in the H-T-plane as observed experimentally is 
shown in Fig. [TJ In addition to the phases shown, a pos- 
sible second transition and spin-liquid phase at H = 
just above T c has been reported recently^ 

Another interesting aspect of MnSi is its sensitivity 
to hydrostatic pressure. With increasing pressure p, 
the magnetic transition temperature T c decreases until 
it vanishes at p — p c rj 15kbaiX£ The transition is 
second order or very weakly first order above a tem- 
perature of approximately 10 K, and strongly first or- 
der at lower temperatures, with a tricritical point sep- 
arating the two regimes. These features have been ex- 
plained as universal properties of quantum ferromagnets 
in an approximation that neglects the helical order at 
longer length scales^ In the paramagnetic region at 
pressures up to approximately 2p c and at temperatures 
T < 10 K, strong non- Fermi-liquid transport properties 
are observed, with the electrical resistivity p c \ displaying 
a temperature dependence p c \ oc T 3 / 2 over almost three 
decades in temperature^ The origin of this behavior has 
recently been proposed to be a combination of columnar 
spin textures and very weak quenched disorder.— In a 
smaller region, above p c and below a pressure-dependent 
temperature To, short-range (spin- liquid) helical order 
has been observe d 17 ' 18 and various explanations in terms 
of analogs of blue phases in liquid crystals have been 
proposed*^ - — Recently, an anomalous Hall effect has 
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FIG. 2: Schematic phase diagram of MnSi in the T-p plane. 
The tricritical point (TCP) separates a line of second order 
transitions (solid line) from a line of first-order transitions 
(dashed line). The inset shows the resistivity data from Ref. 
[la in the non-Fermi-liquid (NFL) region. The boundary of 
the NFL region (dotted line) is not sharp. SL1, SL2, and SL3 
refer to the possible spin-liquid phases or regions reported in 
Refs. fl7l and [Isl. l23l. andfiH. respectively. 




1 st order 



FIG. 3: Schematic phase diagram of MnSi in the space 
spanned by T, p, and H showing the tricritical wings and 
quantum critical points (QCP). 



2. Properties of the helical phase 

If one ignores crystal- field effects, i.e., in an isotropic 
model, the long-range helical order in the helical phase 
(H = 0, T < T c , p < p c ) breaks the translational sym- 
metry. This spontaneous breaking of a continuous sym- 
metry leads to a Goldstone mode, the helimagnon^ The 
dispersion relation is anisotropic, and given by 



whmO) = \J c \\ k f\ +c ± k 4 L . (1.1) 

Here the wave vector k = (ku,k±) has been decomposed 
into components parallel and perpendicular, respectively, 
to the helical pitch vector q, and cm and cj_ are elastic 
constants. The absence of a term proportional to k^_ un- 
der the square root in Eq. (jl.ll) reflects rotational sym- 
metry and is analogous to the nature of the Goldstone 
mode in either smectic or cholesteric liquid crystals^ 
Crystal-field effects that couple the electrons to the un- 
derlying lattice break the rotational invariance and lead 
to a k\ term with a small pref actor, i.e., the dispersion 
relation becomes phonon-likc. Both the single-particle 
relaxation rate 1/r that results from the scattering of 
electronic quasi-particles by helimagnons and the ther- 
mal resistivity pth have a non-Fermi-liquid temperature 
dependence 1/r oc p t h oc T 3 / 2 . In contrast, the transport 
relaxation rate that determines the electrical conductiv- 
ity p c i in a Boltzmann approximation has a temperature 
dependence 1/r^f oc p c \ cx T 5 / 2 . 2 ! The properties of the 
helical phase in the presence of weak quenched disorder 
have also been investigated , 28 i 29 

In the present paper we will discuss the derivation of 
the phase diagram sketched above, as well as some re- 
finements, within the context of a mean-field theory and 
with an emphasis on the H-T plane. We will also deter- 
mine the nature of the Goldstone modes, and their con- 
sequences for electronic properties, in various observed 
or proposed phases with long-range order. In addition 
to the specific heat, the single-particle relaxation rate, 
and the electrical conductivity we will consider the ther- 
mal resistivity pth, whose temperature dependence is the 
same as that of the single-particle relaxation rate. 



B. Model, and energy scales 



been observed at intermediate pressures and below T c 
in the helically ordered phase, and it has been suggested 
that this feature is related to the short-ranged order ob- 
served above p c ^ The phase diagram in the T-p plane 
is shown in Fig. [2j If a magnetic field is applied in the 
vicinity of p c , tricritical wings, i.e., surfaces of first-order 
transitions, emerge from the tricritical point that are be- 
lieved to end in a pair of quantum critical points in the 
T = plane. 15 This feature, which is depicted in Fig. |31 
has been explained theoretically in Ref. wA. 



1. LGW functional 

We will consider a Landau- Ginzburg- Wilson (LGW) 
functional for a three-dimensional order parameter (OP) 
field M = (Mi, M2, M3) whose expectation value is pro- 
portional to the magnetization. We will consider an ac- 
tion that is appropriate for the helical magnets MnSi and 
FeGe, which crystallize in the cubic B20 structure with 
space group P2i3. We will organize the action according 
to the dependence of its various constituents on powers 
of the spin-orbit interaction g so . Within this scheme, we 
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write 

where 
Ah 



A = Ah + A 



DM 



A cf , 



(1.2) 



dx 



t -M\x) + ^{VM{x)f 
d - (V-M(x)) 2 + |(M 2 (x)) 2 



~(VM 2 (x)) -H-M(x) 



(1.3) 



is the action for an isotropic classical Heisenberg ferro- 
magnet in a homogeneous external magnetic field H. 
J v dx denotes a real-space integral over the system 
volume. (VM) 2 stands for £i,j=i Mj d i M j with 
<9i = d/dxi the components of the gradient operator 
V = (81,82,83) = (d x ,d yi d z ). t, a, d, u, and w are 
the parameters of the Landau theory and are related to 
the microscopic energy and length scales (see Sec. II B 21 
below); they are of zeroth order in the spin-orbit coupling 
g so . Equation (|1.3[) contains all analytic terms invariant 
under co-rotations of real space and OP space up to quar- 
tic order in M and bi-quadratic order in M and 
We have added one higher-order term, with coupling con- 
stant w, as an example of a class of terms that can sta- 
bilize unusual phases in helimagnets, although they are 
not of qualitative importance in ferromagnets. We will 
consider this term only in Appendix |A"1 



*4dm = £ / dx M(x) • (V x M(x)) 
2 Jv 



(1.4a) 



is the chiral Dzyaloshinksy-Moriya term that favors a 
nonvanishing curl of the magnetization. The existence 
of this term hinges on the spin-orbit coupling, as well as 
on the system not being invariant with respect to spatial 
inversion (due to the linear dependence on the gradient 
operator). The coupling constant c is linear in g so , and 
on dimensional grounds we have 



c = ak F g s 



(1.4b) 



with the Fermi wave number which serves as the mi- 
croscopic inverse length scale. In this context, this can 
be considered the definition of g so . 

The preceding contributions to the action are all invari- 
ant under either separate rotations, or co-rotations, in 
spin (or magnetization) space and real space. The spin- 
orbit interaction couples the electron spin, and hence the 
magnetization, to the underlying lattice. Therefore, in 
addition to the rotationally invariant terms, any term 
that is invariant under elements of the space group con- 
nected with the crystal lattice is allowed. For the B20 
structure of MnSi and FeGe, the appropriate space group 
is P2i3. To quartic order in M, and bi-quadratic order 



in V and M, the allowed terms in the action are the 
crystal-field terms 



Ac 



r 3 r 



h -(d iMi { X )f + h -^{d i M i+1 {x)f 



(1.5a) 



where M4 = M\. The last term is the usual cubic 
anisotropy that is always present in a magnet on a cubic 
lattice, and 



1 4 

v = u 9so 



(1.5b) 



with \u'\ « u. Of the gradient-squared terms, the first 
one also has cubic symmetry; the second one does not, 
but is invariant under elements of P2i3. On dimensional 
grounds, we have 



b 



a'g 2 a , 
a'i9so, 



(1.5c) 



with I a' I 



a. 



2. Length and energy scales 

The various contributions to the action A, and their 
dependencies on g so , imply a hierarchy of energy scales 
and corresponding wave number or length scales. At ze- 
roth order in <? so , we have the microscopic scale, which is 
represented by the Fermi energy ep and the Fermi wave 
number fcp. Fluctuations renormalize this to the criti- 
cal scale, which is represented by the magnetic ordering 
temperature T c and the corresponding length scale. The 
physics at these scales is described by „4h, Eq. (|1.3[) . 

The chiral DM term is balanced by the rotationally 
invariant gradient squared term in Eq. (11.31) that makes 
magnetization gradients energetically costly. As a result, 
the relevant gradient or momentum scale is of O(g so ), and 
hence the chiral wave number scale is given by the micro- 
scopic scale times g so . This determines the parameters 
of the helix, in particular the helical pitch wave number 
q oc g so . In MnSi and FeGe, this wave number scale is 
on the order of 100 times smaller than the microscopic 
scale. Adm contains one explicit factor of g so and one 
gradient, and hence its contribution to the free energy is 
of O(<7 2 ). The physics at this scale is described by ^dm 
in conjunction with Ah- 

At fourth order in g so , crystal-field effects come into 
play. They pin the helix, are small compared to the chiral 
energy scale by another factor of g 2 , and are described by 
Act, Eqs. (|1.5[) . Since gradients are effectively of O(g so ), 
see above, the contributions of all three terms in A c f to 
the free energy are of 0(gf o ). 

Finally, the external magnetic field sets a scale that is 
continuously tunable. 
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II. PHASE DIAGRAM 

We now derive the mean-field phase diagram for sys- 
tems described by the action given in Eqs. (jl.2p - (|1.5I) . 
We will use the hierarchy of energy scales explained in 
Sec. lIBI to show how a more and more sophisticated phase 
diagram emerges as one keeps effects of higher and higher 
order in g so . 

To do so, we consider field configurations of the follow- 
ing form: 

M(x) — mo+m| e + cos(g-a;)+m^ e sin(q-a;). (2.1a) 

Here mo is a homogeneous component of the magneti- 
zation, rri^ are amplitudes of Fourier components with 
wave vector q, and e± are two unit vectors that form a 
right-handed dreibein together with q: 



qxe+ 



xq 



(2.1b) 



where q = q/q. The sinusoidal terms in Eq. (|2.1a[) de- 
scribe a helix with pitch vector q. The helix is in general 
elliptically polarized, and it is useful to define a polariza- 
tion parameter 



7r = m l jm\ , 



(2.2) 



Special cases are circular polarization, 7T = 1, and linear 
polarization, tt — or ir — oo. The motivation for the 
ansatz Eqs. (|2.1[) is provided by the fact that it gives the 
functional form of the global minimum of the action A, 
Eqs. (|1.2[) - (|1.5p . if one neglects the crystal- field terms 
Ad] that is, for the action up to 0(g 2 o ), see Ref. [Ill and 
Sec. [OB] below. 



A. 0((?s ): Ferromagnet 

To zeroth order in g so the system is approximated 
by a ferromagnet. According to the action An, Eq. 
(|1.3[) . for H = there is a second-order phase transi- 
tion which in mean-field approximation occurs at t = 0. 
For H there is a crossover at t = from a field- 
polarized paramagnetic state, where the magnetization 
extrapolates to zero for H —> 0, to a field-polarized ferro- 
magnetic state where the magnetization extrapolates to 
mo = -\/—t/uH. The free energy density in mean-field 
approximation, / = A/V, in a zero field is 



In a nonzero field, one has 



(2.3) 



(2.4a) 



where mo is the solution of the mean-field equation of 
state 



This is just a classical Heisenberg model, so it cannot 
explain the tricritical point and the associated tricritical 
wings. The latter emerge within a renormalized mean- 
field theory that takes into account the coupling of the 
magnetization to other electronic degrees of freedom^ 
This leads to a fluctuation-induced first-order transition 
in analogy to the case of the nematic-to-smectic- A transi- 
tion in liquid crystals^ as well as to the tricritical wings 
and the quantum-critical points in an external magnetic 
fields The ferromagnetic approximation suffices for un- 
derstanding the gross features of the phase diagram in 
the T-p plane. 



B. O(ffso) : Helimagnet, conical phase 

To second order in <7 so we need to add the DM term, 
Eq. (|1.4a[) . to the action. It is obvious that this term 
favors a nonzero curl of the magnetization, with the di- 
rection of the curl depending on the sign of c. However, 
the spatial variation of M will be limited by the other 
gradient terms in the action, the (VAf ) 2 term in partic- 
ular. We thus expect a spatial modulation of M on a 
length scale on the order of a/c. It is easy to check that 



the ansatz, Eqs. (|2.ip . with n = 1, i.e., 



mi, 



and 



m 

m 

m? 



qH, 
m H, 

H/{cq-aq 2 ) 
1 
u 



—t + cq — aq — 



H 2 



(cq — aq 2 ) 2 



(2.5a) 
(2.5b) 
(2.5c) 

(2.5d) 



solves the saddle-point equations for the action An + 
.4dm- In order to determine q, we extremize the resulting 
free energy with respect to q. q = c/2a is a solution for 
all values of H . Finally, we need to ascertain that the 
solution is a minimum, which turns out to be the case 
for t < aq 2 and H < aq 2 ^—(t — aq 2 )/u. We thus find 
that the field configuration 



M{x) = mo H+mi e + cos(qH ■ x) + e sin(qH ■ x) 

(2.6a) 

with ei, e 2 , H forming a dreibein, and 



m 
mi 



c/2a, 
H/aq 2 , 



where 



= V-r/u y/i - (H/H c2 y 



r = t — aq , 
H c2 = aq 2 y / -r/u, 



(2.6b) 
(2.6c) 
(2.6d) 

(2.7a) 
(2.7b) 



t mo 



H. 



(2.4b) 



minimizes the free energy in the parameter range r < 0, 
H < H C 2. We will refer to this state as the aligned coni- 
cal state (ACS), to distinguish it from the perpendicular 
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conical state discussed in Appendix [A] The ACS is ac- 
tually a global minimum, as can be seen by writing the 
action as a sum of positive semi-definite terms that are 
individually minimized by this stated The mean-field 
free energy density in that range is 

-1 fr 2 H 2 \ 

Equations (|2 .6[) describe the helical phase for H = and 
the conical phase for < H < H C 2- Comparing Eqs. (|2.3p 
and ()2.8[) we see that the helical transition preempts the 
ferromagnetic one. 

For H — > H C 2 from below the helical component of the 
magnetization vanishes, and the free energy, Eq. (|2.8|) . 
approaches that of the ferromagnet, Eqs. (|2.4[) . For H > 
H C 2 the equation of state and the free energy for the DM 
action An + ^4dm are the same as for a ferromagnet with 
action An- 

These considerations account for the structure of the 
phase diagram shown in Fig. [T] except for the field H c i 
and the A-phase. 

C. 0(gl ): Helimagnet, A-phase 




FIG. 4: Spin configuration as given by Eq. (|2.9[) with the 
q vectors in the x-y plane, q = 0.0133A -1 , mo = 0.0146, 
and mi = 0.0323. The arrows represent projections of the 
spins into the plane. The dark and light regions denote spin 
directions antiparallel and parallel to z, respectively. 



Recent neutron-scattering experiments by Miihlbauer 
et ali^i are consistent with the notion that the A-phase 
is characterized by spin textures that form line defects in 
the direction of the magnetic field, with the lines forming 
a hexagonal lattice, the A-crystal, in the plane perpendic- 
ular to the field. This experimental observation led the 
authors of Ref. [ll| to suggest a skyrmion state consisting 
of three co-planar helices as underlying the A-phase: 

M(x) = m H 



-mi 



E 

i=l 



[e^ cos(g( 4 ) • as) + e [ l ] sin(gW ■ x)] , (2.9) 



where the pitch vectors q^' all have the same modulus 
q, are perpendicular to H, the angle between adjacent 
pitch vectors is 2ir/3, and , e_ , and form a right- 
handed dreibein for each value of i. The values of q, mo, 
and mi are obtained by minimizing the action. 

The free energy difference between this state and the 
conical one has a minimum at H ss 0AH C 2, but it is 
still positive even at the minimum. However, Ref. [ll| 
found that taking into account Gaussian fluctuations 
stabilizes the state with respect to the conical one. It 
should be noted that Eq. (|2.9[) is not a solution of the 
saddle-point equations for the action given by Eqs. (II. 2^ - 
(jl.5[) , and therefore cannot be a true local minimum of 
the free energy. Also, the relation between this ansatz 
and what are commonly called skyrmionic spin configu- 
rations, which are solutions of the saddle-point equations 
(see Ref. [2l| and references therein), is not clear. Most 
likely it represents a single-Fourier-component approxi- 
mation to a saddle point and is a Bloch-type description 



of skyrmions as opposed to a Wannier-type description of 
isolated skyrmions in Ref. [2l|. In any case, it describes a 
hexagonal array of line defects with the spin antiparallel 
to the magnetic field at the defect centers, and parallel 
at points on the cell boundaries, see Fig. HI in qualitative 
agreement with the neutron scattering data. 

Earlier, Grigoriev et al.— had proposed a single-helix 
state with the pitch vector oriented perpendicular to the 
external field. Although current experimental evidence 
favors a skyrmion state as a candidate for the A-phase, 
it is still of interest to discuss such a perpendicular conical 
state (PCS), since it might be a viable candidate for the 
ground state in some other part of the phase diagram of 
helimagnets. We therefore briefly discuss the PCS and 
its properties in Appendix [AJ 



D. O(#so): Helimagnet, NFL region 

The skyrmion lattice approximately described by Eq. 
(|2.9[) can melt, which will lead to a skyrmion liquid. In 
such a state the line defects shown in Fig. @] still exist, 
but they no longer form a lattice. Rather, their fluctu- 
ations in the plane perpendicular to the line, which are 
illustrated in Fig. [5j have become so large that the long- 
range order is destroyed. Such a state has recently been 
proposed to represent the NFL region shown in Fig. [2j 16 
Although it is not an ordered phase, such a state has 
much in common with the A-phase and we will discuss it 
in the context of the ordered phases. 
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FIG. 5: Schematic rendering of the hexagonal lattice of 
skyrmion lines shown in Fig. [4] and of fluctuations about this 
state, a is the lattice constant, and u(x) is a displacement 
vector. 

E. 0(gt ): Crystal-field effects 

To fourth order in g so , we need to take into account 
the crystal- field terms shown in Eq. (|1.5a|) . This makes 
the saddle-point equations very complicated, and no ex- 
act solution is known. We therefore take a variational 
approach by inserting Eq. (|2.1al) into the action and min- 
imizing with respect to the parameters of the ansatz. Of 
all the members of the class of functions represented by 
Eqs. (|2.1[) this will yield the one with the lowest free en- 
ergy 

By writing Eq. (|1.5al) we have fixed the coordinate 



system by choosing the crystallographic axes to be the 
x,y,z = 1,2,3 axes. We thus are no longer free to choose 
the direction of H, mo, or q. For simplicity, we will con- 
sider only the case of a magnetic field along the z-axis: 
H = (0,0, if), q = q/q we parameterize in terms of 
angles i? and tp as follows, 

q = (sinz? cosp, sin$ siny>, costf) = (/3i, /3 2 , Pa), 

(2.10a) 

with Yli=iPi — 1- This leaves one free parameter for 
e_|_, namely, an azimuthal angle tp e : 

e+ = (cos 1? cos tp, sin ip e — sin tp cos tp e , 
cos $ sin ip> sin tp e + cos tp cos tp e , 
-sintfsin(^ e ). (2.10b) 

This uniquely determines e_ = q x e + . Finally, mo in 
general needs to be decomposed into components parallel 
and perpendicular, respectively, to q£l However, while a 
perpendicular component can lead to a slightly lower free 
energy, it has no qualitative effects on the structure of the 
phase diagram, and we therefore restrict our ansatz to 

m = mo q. (2.10c) 
We further assume that the system is sufficiently close 
to a second order or weakly first order phase transition 
that one can neglect the last term in Eq. (jl.5al) . With 
these approximations, the free energy does not depend on 
the angle tp e and is completely parameterized in terms of 
six parameters, namely: two amplitudes, mo and mi = 
[(mf) 2 + {m^) 2 ] 1 / 2 /^/2, the polarization parameter n, 
the modulus q, and the two direction angles and tp. 
We find 



2\2 



: 5tm + -(m + mi ) + - (m + m x ) 
- m\ q 2 B B {4, <p) + j ml q 2 B ls (§, tp) 



5t + aq 2 — cq + — cq(Sir)^ 



j ml q 2 Sn B a ($, tp) - ^- m 2 q 2 Sir B la (tf, tp) 



Hmo cos-d 

0(g! o 



(2.11a) 



Here we have defined r = t — St, with St arbitrary at 
this point. We also have made use of the fact that we 
know, from Sec. Ill Bl that the physical solution has the 
property ir = l+<57r with Sir = 0(g 2 o ), and have expanded 
in powers of Sir. The angle-dependent functions in Eq. 
(|2.11ap are 



B S ^, V ) = B + (#,tp)±B-(#,tp), 
Bu,a(0,¥>) = Bf(0,<p)±Bi(0,<p), (2.11b) 



where 



2 sin 2 $ sin 2 tp cos 2 tp, 
sin 2 1? cos 2 z? (1 



• 4 

sin tp - 



B+{§,tp) 
B~ \§,tp) 
B+(0,<p) 

f?r (1?, tp) = sin 4 1? sin 2 tp + cos ■& cos 2 tp 
+ sin 2 $ cos 2 1? sin 2 tp cos 2 tp. 



cos tp), 



sin 2 1? cos 4 tp + cos 2 i9 sin 2 tp, 



(2.11c) 



We now need to minimize the free energy with respect 
to the six parameters. We will first consider the case 
H = to understand the pinning of the helix by the 
crystal-field terms, and then determined the effects of a 
magnetic field. Furthermore, in order to keep the discus- 
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sion transparent we will initially restrict ourselves to an 
ansatz with a circular polarization, Sir — 0. This suffices 
to understand the existence of the critical field H c x- We 
will then generalize the ansatz to allow for a non-circular 
polarization and show that this leads to a splitting of 
the transition at H c \, with a first order transition from a 
circularly polarized helix to an elliptically polarized one 
at a critical field H c i < H c i preceding the alignment 
transition at H c i. 



where 



with 



St = 



t-St, 



74 



(6 + W 

These results are valid for H = and b < bi/2. 



(2.13b) 



(2.13c) 



1. Pinning of the helix 



2. The alignment transition, and the critical field H c i 



We consider the system at H — and initially restrict 
our ansatz to the case of circular polarization, Sn = 0. 
The remaining angular dependence in the free energy, 
Eq. (|2.11a[) . is contained in the functions B s and B\ s . 
Minimizing with respect to ip we find that there are two 
cases. 

Case (1): tp — 0, tt/2, tt, 3it/2. This implies q = 
(0, ±sinz9,cos?J) or (± sinz?, 0, cosi?), and B s = 
i sin 2 2??, Bis — sin 2 i9 + cos 4 i9. Minimizing with 
respect to leads to two subcases: 

Case (la): 2b < b\. The free energy is minimized 
by i9 = 7r/4, 37r/4, 5ir/4, 7tt/A, which implies 
q = (1, 0, l)/\/2 or equivalent. 

Case (lb): 2b > b±. The free energy is mini- 
mized by $ = 0, 7r/2, 7r, 37r/2, which implies 
q = (1, 0, 0) or equivalent. 

Case (2): <p = 7r/4, 37r/4, 57r/4, 7ir/4. This implies q = 
(sin^/v^, sin^/v^, cosd) or equivalent, and B s = 
1 - sin 4 - cos 4 <&, B ls = 1 - sin 2 d + 3 sin 4 
Minimizing with respect to ?? yields 

Case (2a): 2b < bi. The free energy is mini- 
mized by = ± arcsin ^/2/3, which implies 
q = (1, 1, l)/\/3 or equivalent. 

Case (^J: 2b > b x . The free energy is minimized 
by 1} = 0, 7r, which implies q = (0, 0, ±1). 

By comparing the resulting free energies for these cases 
we see that Case (lb) provides the minimum for b > &i/2, 
whereas Case (2a) provides the minimum for b < 6 1 J2 . 
This is a generalization of the result obtained in Ref. l33l 
which considered a model with 61 — 0. 

In MnSi, the pinning is observed to be in the (1, 1, 1)- 
directions, which implies b < b\/2, and we will focus on 
this case from now on. Minimizing the free energy with 
respect to q, we find 



Q = 



c/2 



a+(6 + 6i)/3' 



(2.12) 



which generalizes Eq. (I2.6b[) . Choosing St = cq/2 and 
minimizing with respect to mi we finally have 



For H > we expect the pitch vector to move away 
from (1, 1, 1) towards (0, 0, 1). The calculation proceeds 
as for H — 0, except that now the minimization with 
respect to ■& yields an iJ-dependent result. For Case (2a) 
we find 

a o Jl liH>H cl 

\\{l + 2H 2 /H cl 2 ) UH<H cl , 

(2.14a) 

where (remember r < and b < b±/2 in the ordered 
phase for Case (2a)) 



H, 



cl 



r{b- b 1 /2)(Stf/ua. (2.14b) 



/ = -r 2 /4 W , 



(2.13a) 



The helical pitch vector q thus moves continuously along 
the shortest path from its initial value, (1,1, l)/v3 at 
H = 0, to (0,0,1) at H = H c x, and remains in that 
position for H > H c i. There thus is a second order 
transition at H — H c i— that we refer to as the alignment 
transition. An inspection shows that Case (lb) has a 
larger free energy for all H < H c \. 

For reference in the next subsection we mention that 
if one expands the free energy for small values of 1?, and 
looks for the instability of the solution with § = 0, one 
finds that the latter occurs at H = H cl , as expected. 



3. The polarization transition, and the critical field H c / 

The circular polarization ansatz we have used so far 
explains the two critical fields H c \ and H C 2 observed in 
MnSi. However, the solution obtained in this way misses 
a qualitative feature, as was first pointed out in Ref. [34] on 
symmetry grounds. Since in general b ^ b\, the crystal- 
field contribution to the action, Eq. (|1.5a[) . is not invari- 
ant under x y. As a result, there is no reason for the 
x and y-components of q to become nonzero at the same 
value of H as H is lowered from above, yet the solution 
constructed in the previous subsection has this property. 
Clearly, this is a result of the fact that our ansatz with 
circular polarization possesses cubic symmetry, while the 
action does not. In general, one therefore expects two 
separate transitions in the vicinity of H c \] one where the 
x-component of q becomes nonzero, and a separate one 
where the y-component becomes nonzero. We now show 
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that this expectation is borne out if we allow for a non- 
circular polarization of the helix, which breaks the cubic 
symmetry of the ansatz. 

Consider the full Eq. (|2.11al) . allowing for Sir ^ 0. Min- 
imizing with respect to ip we see that there are two dis- 
tinct cases. 

Case (1): tp — 0,7r /2,7r , 37r /2 and Sn arbitrary. Minimiz- 
ing with respect to Sir yields 



f il 



2c 



(2.15a) 



That is, the polarization is in general elliptical. We 
have B + = 0, B~ = 2 sin 2 $ cos 2 z?, which leads to 



Bs,a = ±^ sin 2 2^9, 
S ls = sin 2 tf + cos 4 tf. 



(2.15b) 
(2.15c) 



Considering B\ & we find two subcases. The first 
one is 

Case (l)(i): p — 0, n, which implies 

q = (0,± situ?, cos 0), (2.16a) 

and 

B u = sin 2 - cos 4 0. (2.16b) 
The second one is 

Case (l)(ii): tp — 7r/2,37r/2, which implies 

q = (±sin??,0,costf), (2.17a) 

and 

a = cos 2 d ± sin 4 §. (2.17b) 

The second case is 

Case (2): tp = tt/4, 3tt/4, 5tt/4, 7tt/4 and 

St: = 0. (2.18a) 

That is, the polarization is circular. This implies 

q = [ — sin i?, — sini9,costf ) , (2.18b) 
Vv2 V2 / 



and 



B s = 1 - - sin 4 - cos 4 1?, (2.18c) 
B ljB = 1 - sin 2 + ^ sin 4 i? (2.18d) 



Now first consider the case H = 0. Minimizing with 
respect to $ one finds that, for b < b\/2, Case (2) yields 
the lower free energy^ The physical solution is thus a 
circularly polarized helix pinned in the (1,1,1) directions, 



(2) 
(1) 



J I 



H'H, H 

cl cl 



FIG. 6: Schematic plot of the free energy density as a function 
of H for Case (1) (dashed line) and Case (2) (dotted line), 
respectively. 



and the relaxation of the condition we had imposed in 
Sec. Ill EH does not change anything. For b > &i/2, the 
physical solution is an elliptically polarized helix pinned 
in the (0, 0, 1) directions. 

Next we consider a solution with $ — 0, which we 
expect to be stable for sufficiently large H . It is easy to 
see that Case (1), which takes advantage of the possibility 
of an elliptical polarization, has a free energy that is lower 
than that of Case (2) by a term of 0(g® o ) everywhere in 
the conical phase, i.e., for H < H C 2- We next look for 
the instability of the & = solution at small H, which 
can be found by expanding the action to second order 
in ■& and looking for the zero of the coefficient of the 
quadratic term. As expected, this instability occurs at 
a field H c i = O(g^ ). To leading order in g so , H c i is 
given by Eq. (|2.14bl) . The values of H c i for the two cases 
are different, but the difference is only of O(g^ ) and is 
irrelevant for the following argument. 

We now have the following situation. For H < H c \, 
the free energy of Case (1) is lower by a term of 0(g® o ). 
However, we know that at H = the free energy of Case 
(2) is lower by a term of O(g 4 ). The two solutions cross 
at a field H c \ given by 



1 



1 



3bl 



2 y Aa\b-h/2\ 



1 - 0{g s 



(2.19) 



see Fig. [B] At this value of H, the state will change 
discontinuously from an elliptically polarized helix with 
a pitch vector given by either Eq. (|2. 16a|) or (I2.17a|) to a 
circularly polarized one with a pitch vector given by Eq. 
(l2~T8bl) . 

We now have the following progression of phases and 
phase transitions as the magnetic field is lowered from a 
value greater than H C 2~. 

H > H c2 : Field-polarized state, no helix. 
H = H C 2 oc g 2 Q : Second-order transition to a conical 
state with an elliptically polarized helix, q — (0,0, 1). 
H = H c i cx <7g : Second-order transition to a coni- 
cal state with elliptical polarization as above, but q = 



10 



H 
H 



c2 



H 
H 



field-polarized FM . 






. field- 

polarized 


\J PM 


conical \ 




I a\ 




pinned helical ""^^^ 





4- Pinning of the A-crystal and the NFL region 

Experimentally, the A-crystal is found to be pinned 
very weakly, and theoretical considerations conclude that 
the pinning potential is only of 0(g^ o )J^ This is a conse- 
quence of the six-fold rotation symmetry of the A-crystal. 
Similarly one expects the (average) direction of the fluc- 
tuating defect lines in the NFL region to be pinned only 
very weakly. 



F. Beyond classical mean-field theory 



FIG. 7: Phase diagram of MnSi in the H-T plane showing the 
predicted polarization transition at H c \ (first order, dashed 
line) in addition to the second-order transitions at H c i and 
H C 2- The nature of the transition from the conical phase 
to the A-phase, or the A-phase to the paramagnetic phase, 
has not been investigated so far. The dotted line denotes 
the crossover from a field-polarized paramagnet to a field- 
polarized ferromagnet. 

(0, sin$, cost?). $ increases from zero with decreasing H. 

H = H c i = H c i[l — O{g so )]: First-order transition 
to a conical state with circular polarization and q = 
(-^ sinz?,-^ sinz?, cos ■&). d increases from its value at 

H c \ with decreasing H . 

H = 0: System reaches helical state with circular polar- 
ization and q = (1, 1, 1)/^. 

The phase diagram is thus predicted to have the struc- 
ture shown in Fig. [JJ with the second order alignment 
transition at H c \ followed by a first-order polarization 
(and re-alignment) transition at H c \ . The latter has so 
far not been observed experimentally. It is an explicit 
realization of the type of transition first predicted by 
Walker— on symmetry grounds. It needs to be stressed 
that the detailed features of this transition are restricted 
by our ansatz, Eqs. (|2.ip ; the states obtained are not true 
minima of the action. However, the basic physical idea, 
which is expected to be realized in the true ground state 
as well, is as follows. In the conical phase, where the he- 
lical pitch vector is aligned with the magnetic field, the 
system can take advantage of the lack of cubic symme- 
try of the action, Eq. (I1.5a[) . by forming a helix with a 
non-circular polarization. This leads to an energy gain, 
compared to a circularly polarized state, of order g® . 
However, it forces the pitch vector into either the y = 
or x — plane, i.e., Case (1) above. At low fields, these 
states have a free energy that is larger by a term of O(<? 4 ) 
than the states with the pitch vector along a cubic diag- 
onal, which requires a circular polarization. The compe- 
tition between these two effects leads to the first order 
transition at H c ±, where both the polarization and the 
orientation of the pitch vector change discontinuously. 



We have treated the phase transitions that we have 
discussed in this section within classical mean-field the- 
ory, and the question arises what changes will result from 
the consideration of fluctuations, classical or quantum. 
One example of the effects of fluctuations is the tricrit- 
ical point and the associated tricritical wings shown in 
Fig. [3J which result from quantum fluctuations that 
can be taken into account within a rcnormalized mean- 
field theor y 14 ' 31 Elsewhere in the phase diagram, fluc- 
tuations are also of qualitative importance. Consider, 
for instance, the transition from the conical phase to the 
field-polarized phase at the critical field H C 2, see Fig. 
[T] The transition is characterized by the vanishing of 
the one-dimensional order parameter mi, the amplitude 
of the helix, see Eqs. (|2.6p . The action is invariant un- 
der mi — > —mi, and we therefore expect this transition 
at nonzero temperature to be in the universality class 
of the classical Ising model. 3 - At zero temperature the 
statics and the dynamics will couple ) 37 ' 38 and one ex- 
pects the quantum phase transition to be described by 
Hertz's model^ 8 . with a scalar order parameter. 



III. GOLDSTONE MODES 

Goldstone's theorem states that if a continuous sym- 
metry described by a Lie group G is spontaneously 
broken, with the remaining subgroup in the broken- 
symmetry phase being H (not to be confused with the 
magnetic field), then in the ordered phase there are n 
soft or massless modes, with n equal to the dimension- 
ality of the coset space All of the ordered states 
discussed in Sec. [IT] break a continuous symmetry, and 
therefore there must be one or more soft modes in the 
ordered phase. 40 The number and functional form of the 
Goldstone modes in the various phases can be determined 
from general arguments. In this section we will do so, and 
we will augment these arguments by explicit calculations 
in some cases where this is feasible. Via a coupling to 
the conduction electrons, the Goldstone modes have in- 
teresting consequences for various observables; this is the 
topic of Sec. [TV] 
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A. O(3so) : Ferromagnons 

To zeroth order in g so the system is a ferromagnet, see 
Sec. Ill Al The relevant symmetry is the rotation group, 
G = 0(3), which in the ordered state is spontaneously 
broken to H = 0(2). Hence there are dim(O(3)/0(2)) = 
2 Goldstone modes, the well-known ferromagnons (FM). 
Their dispersion relation for small wave numbers is 



WFM(fc) = Dk 2 



(3-1) 



with D the spin wave stiffness, which depends on the 
magnetization and vanishes as the magnetization goes to 
zero. 

The easiest way to derive this result is to con- 
sider a nonlinear sigma-model-version of the Heisenberg 
action^ Neglecting fluctuations of the amplitude of the 
magnetization, which can be shown to be massive, one 
parameterizes the order parameter 



M (x) = mo 



TT 1 (x) 

7r 2 (a;) 



(3.2) 



and expands the action to bilinear order in 7ri,2- If we 
neglect the small relativistic term with coupling constant 
d in Eq. (|1.3[) , the resulting quadratic form has two iden- 
tical eigenvalues 



A 



^M 2 



m 



(3.3) 



q = H), so we expect one Goldstone mode. At H = 0, 
one expects the soft fluctuations in the ordered phase 
to be phase fluctuations of the helix, and one might 
naively expect them to be governed by an action S — 
J dx (V0(a;)) 2 , with tj> the phase variable, which would 
lead to a soft eigenvalue proportional to k 2 . However, the 
ordered state is also invariant under rotations of q, which 
can be written as a phase fluctuation with a nonvanish- 
ing gradient, so this cannot be the correct answer. The 
lowest-order term allowed by rotational symmetry that 
involves the gradients perpendicular to q is of the form 
(V^it) 2 , with u a generalized phase variable, and this 
leads to an eigenvalue, or inverse order parameter suscep- 
tibility, proportional to the term under the square root in 
Eq. (jl.ip . 42 The dynamics again require additional con- 
siderations, which lead to a resonance frequency that is 
proportional to the square root of the inverse susceptibil- 
ity, unlike the ferromagnetic case, and this is expressed 
in Eq. (jl.ip 4^ These results were first obtained by means 
of an explicit calculation for both classical and quantum 
helimagnets in Ref. |25j. 

A magnetic field breaks the rotational invariance, so a 
fe^-term will be present in the dispersion relation. The 
prefactor is expected to be analytic in H, and we thus 
expect for the dispersion relation of this "conimagnon" , 
the Goldstone model of the ACS, 



wacs(*0 = Jc\\ kf, + 



J U2 

-± K ± 



(3.4) 



with c' ± cx H 2 . 



For H = 0, one has A(fc — > 0) — > 0, which reflects 
the two Goldstone modes. Physically, the homogeneous 
transverse magnetic susceptibility diverges. This is the 
static manifestation of the spontaneously broken contin- 
uous symmetry. Determining the dynamics requires an 
additional step. One either needs to solve an appropriate 
Langevin equation within a classical context^* or treat 
the problem quantum mechanically . 38 ^ Either way one 
obtains Eq. (I3.1[) with D oc mo. 



B. O(<?so) : Helimagnons 

When the spin-orbit coupling is taken into account, we 
have the various phases involving helical spin textures 
discussed in Sec. Ill Al For all of these phases the rele- 
vant symmetry is the translational one. Let T be the Lie 
group of one-dimensional translations. Then the action 
is invariant under G = T<g>T<S)T = T 3 . 



1. Symmetry arguments 

Consider the aligned conical state (ACS) discussed in 
Sec. Ill Bl from which the unpinned helical phase emerges 
as a special case at H = 0. T 3 is broken down to 
T 2 (translations in the two directions perpendicular to 



2. Model calculation 

For the action to 0(g 2 o ), Eqs. (|2.6I I2.T|) constitute an 
exact saddle-point solution, and we can perform an ex- 
plicit calculation of the Goldstone mode. A complete 
parameterization of fluctuations about the saddle point 
can be written 



M(x) — (mo + 5mo(x)) 



+ 



mi + Smi(x) 



VT^Hx) 

cos(qz + ip (x)) 
sin(qz + ip (x)) 
tp(x) 



(3.5) 



The first term is the nonlinear sigma model for the ho- 
mogeneous magnetization from Scc. lIII K\ and the second 
one parameterizes fluctuations of the helix in terms of an 
amplitude modulation, a phase (po, and a third compo- 
nent ip. The amplitude fluctuations one expects to be 
massive, and an explicit calculation confirms this, so we 
drop (5mo and 8m\. The field ip(x) is conveniently writ- 
ten 

i{)(x) = ip+(x) cos(q ■ x) + (fl-{x) sin(q • x), (3.6) 
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which ensure that ip± and (po at zero wave number both 
correspond to M at wave number q. Double counting is 
avoided by restricting the theory to wave numbers small 
compared to q. If we use Eq. (|3.5p in the action to 0(g 2 o ), 
Eqs. (|1 .311 .4a|) . and expand to bilinear order in the fluc- 
tuations (0i, 02, </>3, </>4, 4>b) = (^o^+^-^i,^), we ob- 
tain a Gaussian action 



a 2 q 4 
2uV 



(3.7) 



k i=l 



A sketch of the derivation, and the explicit form of the 
matrix 7, are given in Appendix [B] Of the five eigenval- 
ues, one goes to zero for k —> 0, in agreement with the 
expectation from the symmetry arguments given above. 
It takes the form 



Ai = ak z 



5k* 



with k = k/q, and coefficients 



a 



,i,2 



s = - 



hf + ml) 

(1 + mf) 3 — TOq(1 + m\) + 2m\m\ 



(1 



t>2\3 



rn 



0) 



(3.8a) 

(3.8b) 
(3.8c) 

(3.8d) 



Here we have defined to 2 , 1 = um^/aq 2 . This result 
is consistent with Eq. (|3.4|) . For H = 0, which implies 
too = 0, it reduces to the helimagnon result of Ref.H^. In 
addition, there are four massive eigenvalues that appear 
in pairs. At zero wave number, they are 



A 2 = 
A 3 



A 3 = 
= A 4 



A^ 
= A„ 



rh\[l - 
--m 2 [l 



0(mt)], 
hO(mt). 



(3.9a) 
(3.9b) 



We recognize X v as representing the massive helimagnon 
modes, 25 modified by the the presence of Too, and A w 
as the massive (due to the presence of a magnetic field) 
ferromagnons, Eq. p. 31) . modified by the presence of m\. 



C. O(<?so) : Skyrmionic Goldstone modes in the 
A-phase and the NFL region 

The helical states that have been proposed as candi- 
dates for the A-phase and were discussed in Sec. Ill El 
are not saddle-point solutions of the action, which pre- 
cludes a model calculation of the Goldstone modes re- 
sulting from this type of order. However, assuming that 
the order is stabilized by some mechanism, the functional 
form of the soft modes can still be determined by sym- 
metry arguments analogous to those put forward in the 
previous subsection. 

Consider a skyrmion lattice state. The state described 
by Eq. (|2.9|) is invariant only under translations in one 
direction, viz., the direction perpendicular to the plane 



of helices. The same is true for any state that is charac- 
terized by columnar order, so this property does not de- 
pend on the precise nature of the skyrmions. Any such 
state will thus have dim(T 3 /T) = 2 Goldstone modes. 
This was to be expected: Since the skyrmions form a 
two-dimensional lattice, there should be two generalized 
phonon modes, namely, a compression mode and a shear 
mode. In zero magnetic field, the energy would still be 
invariant under global rotations of the skyrmion lattice. 
Hence, the soft eigenvalue can have no k 2 contribution. 
For H 7^ this is no longer true, and we thus expect 



A = ak\ + (3 k 4 + ^k 2 z , 
for the soft eigenvalue, with 7 cx H 2 , and 



^SL(fe) = V c_l k\ + d z k 2 



(3.10) 



(3.11) 



for the dispersion relation, with c±,c z — 0(1) and c' z — 
0(H 2 ). 

If the NFL region can be interpreted as a molten A- 
crystal, see Sec. Ill Dl then one of the two Goldstone 
modes, the compression mode, will persist as long as their 
is columnar structure. This is analogous to the fact that 
longitudinal phonons, or ordinary sound, exist in a liquid, 
whereas in a crystal one has transverse phonons or shear 
modes in addition. The NFL region extends to H = 0, 
where the dispersion relation of the compressional Gold- 
stone mode is given by Eq. (|3 . 11 1) with c' z = to first 
order in g so . 



D- 0(gt ): Effects of the crystal-field terms 

1. Symmetry arguments 

Now consider the crystal-field terms in the action that 
first appear at O(g 4 ), Eq. (|1.5a|) . For simplicity, let us 
consider the term with coupling constant b. It breaks ro- 
tational invariance, which invalidates the argument that 
leads to the absence of a fc^-term in the soft mode reso- 
nance frequency. The system must remain stable regard- 
less of the sign of b, and we thus expect for the dispersion 
relation of the helimagnons in the pinned helical state 



WHM(fc) 



C||fc 2 ! 



■d{k\ 



-c ± k\, 



(3.12) 



with c" oc which replaces Eq. (jl.lj) . For the more 
general model given by Eqs. (|1.2[) - (|1.5p . bi, v, and H 
will also contribute to the elastic constant c" . 



2. Model calculation 

We now check this expectation by means of an ex- 
plicit calculation. For the model with only the first of 
the crystal-field terms present, we have an exact saddle- 
point solution, namely, 

M sp (x) — toi [e+ cos(<7 • x) + e sin(q ■ x)] , (3.13a) 
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with e_|_, e_, and q a dreibein, 



4- A-phase, and NFL region 



mi 



= \J-rju, 



where 



and 



r = t- cq/2, 

(1,1,1)/V3 iffe<0 
(1,0,0) if 6 > 0, 



<l 



c/2{a + b/'i) 
c/2a 



if b < 
if 6 > 0. 



(3.13b) 



(3.13c) 



(3.13d) 



(3.13e) 



The parameterization of fluctuations about this state is 
given by the second term in Eq. (|3.5I 



M(x) = (mi + 5mi) [e + cos(q • x + tfo(x)) 

+e_ sin(q • x + (po(x)) +qip(x)] , (3.14) 



with ip(x) given by Eq. (|3.6p . We again drop the mas- 
sive amplitude fluctuations and expand the action to 
quadratic order in the phase fluctuations. The Gaussian 
action is of the form 



^E E 

k a=0,± 



<p a (k)T a p(k)(pp(k). (3.15) 



The explicit form of the matrix T depends on the sign 
of 6; it is given explicitly in Appendix [C] An inspection 
of the eigenvalues shows that in either case there is one 
eigenvalue A that vanishes as k — > and hence represents 
the Goldstone mode, as expected from the symmetry ar- 
gument given above. To order k 2 ^ and k^_ we find 



' akl + \bk\ + +\{a + b) k\/q 2 



X 



if b > 0, 



6/3) k 2 



\b\kl + \(a + b + 0(gt))ky q * 
if b < 0. 
(3.16) 

In the case b < we have neglected terms proportional 
to b 2 — 0(.9g O ) in the pref actor of k\/q 2 . This result 
agrees with the functional form obtained by symmetry 
arguments alone, Eq. (|3.1I) . and for b = it correctly 
reduces to the result for the isotropic model, Ref. and 
Eq. CO]). 



3. Generalized helimagnons 

We can now summarize the results for the single-helix 
phases discussed above as follows. In the helical and 
conical phases, including pinning effects, there is a single 
Goldstone mode with a resonance frequency 



Wo(fe) 



C|| k 2 \ + c± k\ + c± kj_ , 



(3.17) 



In the A-phase, pinning effects are weaker than in the 
helical phase due to the hexagonal nature of the skyrmion 
lattice, see Sec. Ill E 41 above. When this weak pinning is 
taken in to account, the Goldstone mode is thus given 
by Eq. (j3~TTj) with c! z = 0{g 6 so ,H 2 ). While in MnSi the 
A-phase is observed only in an external magnetic field, 
there is no intrinsic reason why in some other system it 
could not be stable in a zero field. The Goldstone modes 
in such a system would be given by Eq. ([3~TT|) with an 
extremely small c! z . By the same argument, in the NFL 
region at H = we expect a Goldstone mode given by 
Eq. dSHJ with c' z = 0(gl ). 



E. Summary of Goldstone modes, and temperature 
regimes 

To summarize, we have found that the (single) Gold- 
stone mode in the pinned helical and conical phases is 
given by Eq. (|3.17l) . In the pinned helical phase the par- 
allel direction is determined by the crystal-field effects 
that pin the helix; in the conical phase, it is the direc- 
tion of the magnetic field (which we have chosen to be the 
z-direction for all of our considerations). The elastic con- 
stant c± is small compared to the other elastic constants. 
In the pinned helical phase it is of O(g 2 ), and in the con- 
ical phase it is of 0(H 2 ). By contrast, C|| and c± are of 
O(g so ) and O(g® ), and of O(H ), respectively. In the A- 
phase, there are two Goldstone modes whose dispersion 
relation is given by Eq. (|3 . 11 1) . The elastic constant c' z is 
small of 0(H 2 ) and O(g® ) compared to the others. In 
the NFL region, the single Goldstone mode is also given 
by Eq. (j3~Tl|) . with c' z = O(g s 6 ). Finally, in the perpen- 
dicular conical state discussed in Appendix IA1 there is a 
single Goldstone mode that in the absence of crystal-field 
effects is given by Eq. (|A3|) . If pinning by crystal-field 
terms is taken into account there is an additional term 
Cyk 2 under the square root, with c y — O(g 2 ). All of 
these results are also summarized in Table U in Sec. |Vj 

The structures of the various Goldstone modes, and the 
fact that the various elastic coefficients have very differ- 
ent magnitudes, leads to the formation of different tem- 
perature regimes that are dominated by different physics. 
We now explain this using the generalized helimagnons, 
Eq. (|3.17[) . as an example; the argument for the other 
cases works analogously. As far as the coupling of the 
magnetic Goldstone mode to the electronic degrees of 
freedom is concerned, the resonance frequency luq scales 
as the temperature, ujq ~ T. If we scale fc|| with T/^/cjJ 

1/4 



and k± with VT/c/ , we obtain 



u>o(k)=T x /kf ] 



k±, 



(3.18) 



with c± = 0(H , g 2 D ) small compared to C|| and c±. This 
comprises Eqs. (pO| and (|3.12l) . 



where k now denotes the scaled, dimensionless, wave 
number. For T ^> c±/ \/c± the symmetry-breaking fc^ 
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term is negligible, and the Goldstone mode is effectively 
what it would be in a rotationally invariant system. In 
this regime the physics is dominated by universal hy- 
drodynamic effects that are independent of the micro- 
scopic details of the solid and analogous to the corre- 
sponding effects in liquid crystals. In the opposite limit, 
T <C c±/y/c±, the crystal- field effects due to the un- 
derlying ionic lattice, or the external magnetic field, if 
present, dominate and the Goldstone mode has the same 
functional form as (anisotropic) acoustic phonons. Due 
to the smallness of 5j_ the universal hydrodynamic regime 
is sizable, and it is in this region that the most interest- 
ing effects of the magnetic order manifest themselves in 
the electronic properties of the system. This is true es- 
pecially if the Goldstone mode appears in zero magnetic 
field and the pinning is very small, such as in the NFL 
region or a (so far hypothetical) A-phase in zero field. 



The universal hydrodynamic result, C(T) oc T 2 , was 
first derived in Ref. [27]. It is subleading to, but distinct 
from, the Fermi-liquid result C(T) oc T + 0(T 3 InT). 
At asymptotically low temperatures it crosses over to a 
T 3 behavior consistent with the acoustic-phonon-like dis- 
persion relation in cither the pinned helical phase or the 
aligned conical phase at asymptotically small wave num- 
bers. 



2. A-phase 

In the skyrmion-lattice state we find from Eq. p. lip 
in conjunction with Eq. (|4.ip 

C(T) = const, x i 1/4 v (4.3) 

\T°^/cj_c z ' if T>c;/ v /cl. 



IV. EFFECTS OF GOLDSTONE MODES ON 
ELECTRONIC PROPERTIES 

Via a coupling to the conduction electrons, the Gold- 
stone modes derived in the preceding section influence 
the electronic properties of the helical magnet. In this 
section we derive the consequences for the specific heat, 
the single-particle relaxation rate, and the thermal and 
electrical resistivities. In all cases we consider the con- 
tribution of the Goldstone mode in isolation; it comes in 
addition to all other contributions to these observables. 



A. Specific heat 

Any well-defined (i.e., not overdamped) excitation 
with a dispersion relation Li(k) yields a contribution to 
the specific heat C given by 



C(T) 



d_ 

df V 



(4.1) 



Here Ub(x) = l/(exp(x/T) — 1) is the Bose distribution 
function, V is the system volume, and we use units such 
that H = k& = 1. This allows one to determine the con- 
tributions to the specific heat by the various Goldstone 
modes. 



1. Generalized helimagnons 

We first consider the helical and aligned conical phase. 
The dispersion relation is given by Eq. (|3.17[) . Performing 
the integral in Eq. (|4.1j) yields 

C(T) = const. J T3 J^ , 
V ' \T 2 /^I ifT>c x / A /cI 



(4.2) 



B. Relaxation times and transport coefficients 

The temperature dependence of the single-particle re- 
laxation rate, as well as that of the electrical transport re- 
laxation rate or the electrical resistivity, can be obtained 
by using the results of Refs. [28ll29l The former paper 
obtained the following expression for the single-particle 
relaxation rate, 



1 



r(fc) 



= 2 



du 



-oo sinh(u/T) V 



i^V"(p-fc;fe,p;i 



x8(u-u)i(p)). 



(4.4) 



Here the quasiparticle momentum k is taken to be on the 
1-Fermi surface, wi(fc) = 0, with 



w i,2(fc) = ^ 



6m 



6c ± V (6»+q _ 6 



4A 2 



(4.5a) 

which is separated from the 2-Fermi surface, oJzik) = 0, 
by twice the Stoner gap A. Here 



6; — e fe — £f. 



(4.5b) 



with 

e fc = k 2 /2m e 



— — j (k 2 x k 2 y + k 2 y k 2 z + k 2 z k 2 x ) (4.5c) 



an electronic energy-momentum relation consistent with 
a cubic crystal. The dimensionless parameter v is a mea- 
sure of deviations from a nearly-free electron model. V" 
is the spectrum of an effective potential given by 

V{k; Pl ,p 2 ) = iA 2 X (fc)7(fe,p 1 )7(-fc,p 2 ), (4.6a) 
with vertices 



7(fc,p) 



8to A 



tj (k z p 2 ± + 2(fc± -P±)Pz) 

/vp 

+0(k 2 ), (4.6b) 
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and x the soft-mode susceptibility, 



1. Generalized helimagnons 



X(k) 



2N F 3k* Lj$(k) - (my 



(4.6c) 



Here wo(fe) is the resonance frequency, which is equal to 
the frequency given in Eq. (pR)) . ([53T)l . (|3~T2")l . or ([A"3l) . 
depending on the phase under consideration. The elec- 
trical transport relaxation time, which determines the 
electrical resistivity, is effectively given by averaging a 
similar expression over the Fermi surface,— 



N F 



du 



sinh(w/T) V 



-y 



(p-fc) 2 



xV''(p-fc;fc,p;u)*(u-wi(p))$(wi(fc)). 

(4.7) 

Note the additional, compared to Eq. (|4.4[) . factor of 
(p — fc) 2 under the integral in Eq. (I4.7[) . This is char- 
acteristic of the description of electrical transport in a 
Boltzmann approximation and leads to a temperature 
dependence of the electrical resistivity that is different 
from that of the single-particle relaxation ratei 44 ' 45 In 
contrast, in a Boltzmann description of thermal trans- 
port this additional factor is absent, and the temperature 
dependence of the thermal conductivity is given by that 
of the single-particle relaxation rate4£ 

We note that these expressions for l/r(fc) and l/r*j r 
vanish as q — > 0. While they give the leading asymptotic 
temperature dependence for the relaxation rates in he- 
limagnets, they therefore cannot be used to obtain the 
corresponding results in the ferromagnetic limit. The 
asymptotic low-T behavior of the rates in ferromagnets 
is qualitatively different and briefly treated in Appendix 
ID] In the context of helimagnets, for a non-spherical 
Fermi surface (i.e., v ^ in Eq. ()4.5c[) ^1 . and for generic 
wave vectors k, the temperature scaling behavior of the 
single-particle relaxation rate 1/r = l/r(fc), the thermal 
resistivity pth, and the electrical resistivity p c \ can be 
represented schematically by the expressions 



Pel 



Pth 



dp\\ / dp± 



Pi 



P 



II 



11 J sinh(w (p)/T) 
5(u (p) -p± ~P\\) 



dp I, 



dp 2 ± 



wo(p) 

(p'i+p 



2 \2 

IK 



sinh(w (p)/r) 
6(u (p) ~p± -p\\) 
w (p) 



(4.8a) 



(4.8b) 



The resonance frequency always scales as the tempera- 
ture, wo ~ T, and the temperature dependence of the 
relaxation rates thus is determined by how the momen- 
tum components scale with temperature. 



For the helical and conical phases we have, from Eq. 
(f37TT|) . p|| ~ T, and px ~ T and - T 1 / 2 at asymptot- 
ically low and intermediate temperatures, respectively. 
This yields the following temperature dependence for the 
single-particle relaxation rate and the thermal resistivity: 



1 

T 



T 3 ifT<S_ L / % /7iI 
T 3 / 2 if T>S ± /ycI 



- oc p th oc ^ m3/2 :r ^_ . • (4.9) 



The corresponding result for the electrical resistivity is 

(4.10) 



It 5 ifT«5 ± /^cT 

Pcioc< 5/2 2T>c ± / y /cZ 



In a vanishing external field, and in a temperature regime 
where pinning effects are not relevant, we recover the 
T 3 / 2 and T 5 / 2 behavior for 1/r and p c \, respectively, of 
Ref.i^anddi. 



2. A-phase 

For a skyrmion lattice, we obtain by using Eq. (|3.1ip 

(4.11) 



i It 3 if t« djy/c; 

- OC pth OC { ^ 2 



for the single-particle rate, and 



Pol oc 



(4.12) 



T 5 ifT«c' z /Vc- 
for the electrical resistivity. 

C. Systems with quenched disorder 



The preceding results hold for clean systems. In the 
presence of quenched disorder, elastic scattering of the 
conduction electrons leads to profound effects that man- 
ifest themselves in the transport properties. One needs 
to distinguish between the strong-disorder regime, where 
the transport is diffusive, and the weak-disorder regime, 
where it is ballistic. In a Fermi liquid, these two regimes 
are characterized by Tt -C 1 and Tt 1, respectively^ 
In a helical magnet, the weak-disorder or ballistic regime 
is characterized b y 28 ' 29 



V(e F T i) 2 T/A » 1, 



(4.13) 



where r e i is the elastic relaxation time. In this regime 
there is an additional contribution 5(1 /r) to the relax- 
ation rates that is qualitatively the same for both the 
single-particle rate and the electron transport rate^i and 
thus provides the temperature dependence of the correc- 
tions to both the electrical and thermal resistivities. It 
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was shown in Refs. I28U29I that, for temperature scaling 
purposes, this contribution can be represented by 



6(1/ t) oc 6 pel cx 6p t h 



dun-p(u/T) / dp 1 1 / dp 



± w (p) 



6(u-uj (p)). 

(4.14) 



Generalized helimagnons 



From Eq. (|4.14l) we see that at temperatures where pin- 
ning effects are not important, generalized helimagnons 
lead to 6(1/ V) oc T, a result first obtained in Ref. [29|. 
At asymptotically low temperatures, characterized by 
T <C c±/y/c±, one finds a T 2 -behavior. In the pinned 
helical phase the crossover temperature between these 
two regimes is determined by the strength of the crystal- 
field effects; in the conical phase the magnetic field also 
cuts off the universal hydrodynamics T-behavior. 



2. A-phase 



For the A-phase, Eq. (|3.11j) yields 



r T 2 

5(1/7") « < T 3 /2 



if T « c'JV^ 
if T»c' 2 /V^ 



(4.15) 



The pinning effects in the A-phase are very weak due 
to the hexagonal nature of the skyrmion lattice, with c' z 
only of 0(g® o ), see Secs.|njM|and|niD41 The size of the 
asymptotic region is therefore likely to be dominated by 
the independence of c' z . Whether or not the universal 
hydrodynamic T 3 / 2 -behavior is observable in the A-phase 
(there currently are no experimental indications that it 
is) would require a detailed quantitative analysis that 
goes beyond the scope of the current paper. 



3. NFL region 

The preceding result is also of interest with respect 
to the non-Fcrmi-liquid region shown in Fig. [2J which is 
not a phase with long-range order, but where the electri- 
cal conductivity shows a pronounced T 3 / 2 -behaviorJ£ An 
explanation that has recently been proposed*— is as fol- 
lows. The T 3 / 2 -behavior derived above is a consequence 
of the structure of the Goldstone modes due to colum- 
nar fluctuations, Eq. (|3 . 1 1 1) . in conjunction with weak 
quenched disorder. In the A-phase, which displays long- 
range order in the form of a skyrmion lattice, there are 
two such Goldstone modes, see Sec. IIII CI If the NFL re- 
gion can be interpreted as a melted skyrmion lattice, then 
the resulting skyrmion fluid will still have one Goldstone 
mode with the same structure, namely, the compression 
mode mentioned in Sec. IIII CI Weak quenched disorder 



will then still produce a contribution to the electrical re- 
sistivity, as well as to the single-particle relaxation rate, 
that is proportional to T 3 / 2 in a pre-asymptotic region. 
The NFL region is observed to extend to a vanishing 
external magnetic field, so the low-T boundary of the 
universal hydrodynamic region is determined by the pin- 
ning effects, which are very weak, see Sees. HIE 41 and 
HHP 41 The universal hydrodynamic T 3 / 2 behavior is 
therefore expected to extend to very low temperatures. 
A remaining question is the size of the prefactor, which 
in a bare theory is expected to be small due to the long 
length scale set by the skyrmion lattice. The resolution 
proposed in Ref. [l6| is that mode-mode coupling effects 
drastically enhance the magnitude of the effect, in anal- 
ogy to what is believed to happen in the blue phases of 
liquid crystals ! 48-50 



V. SUMMARY, AND CONCLUSION 

In summary, we have given a comprehensive descrip- 
tion of all phases in Dzyaloshinsky-Moriya helical mag- 
nets where long-range order has been observed. These 
include, the pinned helical phase at weak magnetic fields, 
the conical phase at higher magnetic fields, and the A- 
phase at intermediate magnetic fields and temperatures 
close to the critical temperature. We have shown that 
the system goes from the conical phase at high magnetic 
field to the pinned helical phase at low magnetic field 
via two distinct phase transitions; a second-order transi- 
tion where the orientation of the helix changes smoothly, 
followed by a first-order transition where both the orien- 
tation and the polarization of the helix change discontin- 
uously. For the A-phase we have considered the recent 
interpretation, based on the observation of a six-fold sym- 
metry in the neutron scattering signature, as a hexagonal 
lattice of skyrmionic line defects. We have also discussed 
the perpendicular conical state, which had been discussed 
earlier as a possible realization of the A-phase. In addi- 
tion, we have discussed a proposal for the interpretation 
of the non-Fermi-liquid (NFL) region in the disordered 
phase, which does not display long-range magnetic order 
but has many features in common with the A-phase. For 
all of these states we have determined the number and na- 
ture of the Goldstone modes, and the temperature depen- 
dencies of various observables that result from the scat- 
tering of the conduction electrons by these excitations. If 
the state in question can be described by an exact saddle- 
point solution of a model Hamiltonian we have provided 
an explicit calculation; in other cases we have used sym- 
metry arguments to determine the functional form of the 
Goldstone modes. The results are summarized in Ta- 
ble |U The Goldstone modes show unusual anisotropic 
frequency-momentum relations that result from helical 
spin structures in the various states and are reminiscent 
of the Goldstone modes in smectic and cholesteric liquid 
crystals. The specific heat C(T), the single-particle re- 
laxation rate l/r(T) and the thermal resistivity p^(T), 



17 









ninnpfl Viplipfll 

/conical 


A -nVi ?mp 


PCS" 
(H\\z , q\\x) 


NFL 6 




Goldstone 
modes 


Wo(fe) 










\Jc\\k^ + c ± k 2 ± + c±k\ 








# 


1 


2 


1 


1 


clean 


C(T) 


universal hydro- 
dynamic regime 




rp5/2 




j.5/2 


crystal — field 
regime (T-J-OJ 


r 3 


T 3 


y5/2 


1/t(T), 
Pth(T) 


universal hydro- 
dynamic regime 


ji3/2 




y3/2 




crystal — field 
regime (T-MD) 


r 3 




T 2 


Pd(T) 


universal hydro- 
dynamic regime 


iji5 / 2 


T 3 






crystal — field 
regime (T 1 — yO) 






rj~>3 


weak 
disorder 


Spd(T), 
S P tb{T) 


universal hydro — 
dynamic regime 


T 




T 


j.3/2 


crystal — field 
c 

regime 









"This phase, the perpendicular conical state, has not been exper- 
imentally observed so far. 

''The crystal-field effects are much smaller in the NFL region than 
in the ordered phases, see Sees. IIIE4l and llIID4l We results listed 
are valid in the universal hydrodynamic regime, which is expected 
to extend to very low temperatures; the true asymptotic low-T 
behavior is the same as in the A-phase, where the size of the T — > 
region is determined by the external magnetic field. 

c In the presence of quenched disorder this regime does not rep- 
resent the true asymptotic T — y behavior, which is characterized 
by diffusive rather than ballistic electron dynamics, see the table 
caption and Eq. I|4.13t . 

TABLE I: Properties of various ordered phases and the proposed state representing the non-Fermi-liquid region. Listed are 
the number of Goldstone modes and their respective dispersion relations, as well as the temperature dependence of various 
observables. The elastic constants c_l, c' z , and c z in the Goldstone modes for the pinned helical/conical phase, the A-phase, 
and the PCS, respectively, are due to crystal-field effects or an external magnetic field and are small compared to the other 
elastic constants. The universal hydrodynamic regime is the temperature region where crystal-field effects are not important 
and these elastic constants can be neglected. The true asymptotic behavior as T — > is dominated by the crystal-field effects 
and is realized only at very low temperatures. In the presence of weak disorder, the regime dominated by crystal-field effects 
is bounded below as well as above and does not represent the true asymptotic low-temperature regime. See the text for more 
information. 



and the electrical resistivity p c \(T), all display temper- 
ature dependencies that are distinct from, and in some 
cases stronger than, those in a Fermi liquid. The most 
remarkable result is a non-Fermi liquid T 3 / 2 behavior of 
the electrical resistivity in certain temperature regimes 
in the A-phase and in the NFL region, which is proposed 
as an explanation of the observed enigmatic properties of 
the DM magnet MnSi. 
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Appendix A: The perpendicular conical state 

Guided by earlier experimental results that did not 
reveal the six-fold symmetry discovered by Miihlbauer 
et al.fii Grigoriev et alJ£ had proposed a perpendicular 
conical state (PCS) to represent the A-phase, i.e., a spin 
configuration of the form 

M(x) = moH + mi [e + cos(q • x) + e_ sin(q • x)] 

(Al) 

with q perpendicular to H and e + , e_, and q forming 
a dreibein. This state is not a saddle point either, and 
Ref. [ll] found it energetically unfavorable compared to 
the conical and skyrmion lattice states. However, it can 
be stabilized, at least in principle, and we therefore in- 
clude it in our discussion. Ref. [lfl proposed that the PCS 
is stabilized due to a gap in the helimagnon excitation 
spectrum, which in turn had been proposed theoretically 
in Ref. [5lJ. Such a gap is at odds with Goldstone's theo- 
rem as well as with the calculation in Ref. 25- A different 
possibility, which is consistent with symmetry considera- 
tions, is a large negative value of the parameter w in Eq. 
(jl.3|) , which stabilizes the PCS in a parameter regime 
close to where the A-phase is observed, see Fig. [5] While 
such a large value of w is not realistic, this illustrates 
that the PCS can be stabilized by terms in the action 
that are allowed by symmetry, and it therefore is useful 
to determine the properties of this state. 

Let us now consider the Goldstone modes in the PCS, 
using the symmetry arguments explained in Sec. IHI Al 
For definiteness, we assume that the g-vector points in 
the x-direction. The ordered state is still invariant un- 
der translations in the y and z-directions, so we have 
dim(T 3 /T 2 ) = 1 Goldstone mode. With the magnetic 
field H in the z-direction, the state is still invariant un- 
der rotations of q in the x-y-plane. The eigenvalue there- 
fore cannot have a k 2 term. The soft eigenvalue therefore 
must have the form 



This result is different from the one for the skyrmionic 
Goldstone mode in Sec. IIII CI Measurements of the spe- 
cific heat therefore offer a way to distinguish between 
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FIG. 8: Stabilization of the PCS by a large negative value of 
w. Shown are the regions of stability for various phases for 
the action shown in Eq. (|1.3[) with a = 1, d = 0.025, u — 0.1, 
and w = —300. 



these two states short of a direct determination of the 
spin structure. 

Using Eq. (| A3|) in an obvious generalization of Eqs. 
(I4.8|) we find for the single-particle relaxation rate and 
the thermal resistivity 



i t 2 if T« Cz /v^I r ,,s 



and for the electrical resistivity 



JT 3 ifT«c z /VcT 
A^W i{T>>Cz/v ^ 



(A6) 



\ = akl + f3(k y + ki) + 1 kl 



(A2) 



where 7 cx H 2 for small H. The structure of the dynam- 
ics will be the same as in the aligned conical case, and 
we thus have a dispersion relation 



Appendix B: Gaussian fluctuations in the aligned 
conical phase 



wpcs(*0 = y/Z k 2 x + c z k 2 + c ± (k y + k±), (A3) 

with c z oc H 2 and c x and c± constant for H — > 0. Here 
we neglect pinning effects, which lead to a fc 2 -term under 
the square root at asymptotically low wave numbers. 

For the specific heat contribution we find, from Eqs. 
(ICT) and (TM1), 



C(T) = const, x 



V 5 / 2 / (c x Cz) 1 / 2 C 1 / 4 

T 2 /V^I 



if T«c 2 /7SI 

if T»c,/VcL\ 
(A4) 



Here we derive the Gaussian fluctuation action in the 
conical phase, Eq. (|3.7[) . The equation of state for this 
phase, Eqs. (|2.61 12. 7\i , can be written 



t 



aq 2 + u(rriQ 



mf) 



0, 



(Bl) 



with q and mo given by Eqs. (|2.6bl) and (|2.6c[) . respec- 
tively. Inserting Eq. (13.10)) into the action, Eqs. (|1.3I 
ll.4al) , dropping the amplitude fluctuations Stuq and Smi , 
and expanding to bilinear order in the phases yields Eq. 
(I3~7l) with 
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A determination of the eigenvalues yields Eqs. 



Appendix C: Gaussian fluctuations in the pinned helical phase 



Here we derive the Gaussian fluctuation action in the pinned helical phase, Eq. p,15[) . Dropping the amplitude 
fluctuations, and expanding to linear order in the phase fluctuations, Eq. (|3 . 14[) becomes 



M{x) = M sp {x) + SM(x), 



(Cla) 



with 



SM(x) = mi [(— e_ sin(q • x) + e + cos(q • x)) <po{x) + q sin(q • x) cp-(x) + q cos(q • x) <f+(x) + 0(y>o)] . (Clb) 

M sp is an exact saddle point of the action, Eqs. (|1.2[) - (|1.5|) with b\ = v = 0, so the terms linear in SM vanish. We 
now consider terms bilinear in y>o,+.- and first concentrate on the gradient-free terms. Neglecting rapidly fluctuating 
Fourier components proportional to e tnq ' x with n > 2 one finds 



A {2) [5M]/m\ 



dx (po(x) 



t + aq 2 



t + aq 2 + um\ + bq 2 E] < 



ifa{x) 



+ \ [ dX E V 01 ^ 
J a=± 

— — bq 2 j dx E] e^qf (po(x) <p a (x) + (gradient terms) 

■' a=± i=l 



(C2) 



r 



An explicit calculation shows that 

3 



with 



£a?[(s i + ) a +(*_) a ] = i- mi 

i=l 

E # = 



(C3a) 
(C3b) 



f(q) =pf + (31 + /? 3 4 , (C3c) 
where ^x,2,3 are the direction cosines of q. Furthermore, 

3 



E g « (*) 3 = o 



(C3d) 



for a = +, — and for both q = (l,l,l)/\/3 and q = 
(1,0,0). Using the equation of state, Eqs. p,13b[ I3~.13cp . 
we see that the fo-ifo vertex and the ipo-ip a vertices all 
vanish. At zero wave number we thus have one zero eigen- 
value that corresponds to one Goldstone mode, in agree- 
ment with the expectation from Sec. IIHB II 



We next calculate the gradient-squared terms. Geo- 
metric identities similar to those expressed in Eqs. (|C3p 
allow to determine the vertices. One finds a Gaussian 
action of the form given by Eq. (I3.15[) . For b > 0, the 
matrix T is given by 



J 



r(fc) 



,k 2 + \bk\ 



z 2 c ky 

^ 2 ^ ^ x 



% 2 C ky 



7^ ci k ~\~ tt b k~ 



2 c q i ,, , . ... 2 

— i(a + b)qk z 



% 2 c k x 
i(a + b)qk z 
\cq+ \ak 2 + \bk 2 z 



(6>0) 



(C4) 
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For b < 0, the matrix T takes the form 



r 00 (fe) 


= (a + b/3)k 2 , 








r 0+ (fc) 


= r +0 (fe)* = i (\ 




2v6 




r -(fc) 


= r_ (fc)*=tQ 




6 (fc 2 4 

2\/6 


" - 2fc 2 ), 


r++(fe) 


= r__(fe) = icq 


+ l(a + 6/3)fc 2 , 






r+-(fc) 


= r_+(fe)* = i(a- 


h6/3) 9 fc|| 




(6<0) 



Calculating the eigenvalues one finds Eq. (|3. 161) . 



and for the electrical transport rate 



Appendix D: Relaxation rates in ferromagnets 

In this appendix we sketch how the expressions for 
the single-particle and transport relaxation rates given in 
Sec. UVBl change in the ferromagnetic limit. A more thor- 
ough discussion of this topic will be given elsewhere^ 2 - 

There are three major changes compared to Eqs. (j4.4[) 
and (|4.7p that take place in the ferromagnetic limit. 
Firstly, there is no intra- Stoner-band scattering. This is 
because in a ferromagnet, the Goldstone mode is entirely 
transverse with respect to the direction of the magneti- 
zation, whereas in a helimagnet this is not the c&seM- 
As a result, with k in Eq. (|4.4[) on the 1-Fermi sur- 
face, the resonance frequency inside the delta function 
under the integral will be W2(p). Secondly, the dimen- 
sionless vertex 7, Eq. (|4.6b[) . changes to unity. Techni- 
cally, this is seen most easily within the formalism of Ref . 
l28l . Physically, it reflects the fact that for intra- Stoner- 
band scattering the magnetic fluctuations couple to the 
quasi-particle density, and thus physically act akin to a 
chemical potential, only gradients of which contribute to 
scattering. For inter-Stoner-band scattering, on the other 
hand, they physically act akin to an external magnetic 
field, and thus they couple without gradients. Thirdly, 
the form of the susceptibility Xi Eq. (|4.6c[) . changes. It 
now reflects the ferromagnetic magnons, and for power 
counting purposes it can be adequately represented by 



1 



x"(fe, u) = 1 r r ^5{u - U) FM (k)), 



with 



WFM(fc) = const, x Afc 2 /fcj 



(Dla) 



(Dlb) 



the ferromagnon frequency. Putting all of this to- 
gether, and neglecting numerical prefactors, we obtain 
for the single-particle scattering rate due to ferromag- 
netic magnons 



1 



A 1 



air 



1 



r(fe) N F V ^ sinh(w FM (p)/T) 



S(oj 2 (k+p)-w F u(p)), 
(D2a) 



1 

7t~r 



A 1 



p 2 /k 2 F 



N F V2 sinh ( W FM(p)/T) 



x8(uj 2 (k + p) -u FM (p)) ■ (D2b) 

Performing the integrals leads to the following results. 
For the single-particle relaxation rate on the Fermi sur- 
face, or the thermal resistivity, one finds 



if T < A 3 /e 



- OC pth oc 
r 



Tln(Te 2 /A 3 ) if A 3 /e 2 « T « A (D3) 
Tln(e 2 /A 2 ) ifT>A, 



where we have omitted numerical prefactors as well as 
less-leading terms. For the electrical resistivity, the cor- 
responding result is 



T (A 2 e 2 ' 
/m x <; T 2 /A 
T 



-\*/4t ifT «A 3 /e| 

if A 3 /e 2 « T « A (D4) 
if T » A , 



We provide some brief comments regarding these results, 
a more complete discussion will be given elsewhere^ (1) 
The energy scales A 3 / e| and A that lead to the three tem- 
perature regimes shown above emerge within the bare 
(Stoner-level) theory. It is not obvious how these scales 
become renormalized, and one therefore has to be careful 
when making quantitative comparisons with experiment. 
In particular, within the bare theory the regime T>A 
is not realizable since it implies T > T c ; this may change 
within a fully renormalized theory. (2) The second re- 
sult in Eq. (|D4j) . l/r*j r cx T 2 /A, reproduces the result of 
Ueda and Moriya.— We note, however, that this is not 
the true asymptotic low-temperature result; it is valid 
only in a pre-asymptotic temperature window. (3) The 
true asymptotic behavior for T — > is exponential, rather 
than power-law, in T. This is due to the Stoner splitting 
of the conduction electrons; this effect was neglected in 
Ref. [54|- In helimagnets the situation is qualitatively dif- 
ferent, see Sec. HVBl 
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